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ABSTRACT 


The  allowed  spin  configurations  in  the  rare  earth  orthoferrite 

have  been  derived  on  the  basis  of  group  theory  and  correlated  with  those 

actually  observed.  The  susceptibilities  and  antiferromagnetic  resonance 

inodes  of  YFeO^  have  been  calculated.  The  paramagnetic  resonance  spectra 

of  Fe  and  Gd  in  YAIO^  were  measured  and  the  spectra  analyzed  to  give 

3+^ 

spin  Hamiltonian  parameters  for  these  ions.  The  Fe  site  was  found  to 

3+ 

be  nearly  cublc^  but  the  Gd  site  was  badly  distorted.  The  spin 

Hamiltonian  parameters  were  compared  with  crystal  fields  calculated  on 

a  point  ion  model  using  positional  parameters  measured  not  for  YAlO^  but 

for  GdFeO^.  Poor  agreement  was  found.  Further  evidence  that  the 

crystal  fields  are  substantially  different  in  the  ortho-alumlnates  from 

3+ 

the  orthoferrites  were  obtained  in  the  optical  spectrum  of  Yb  in  YAlO^ 

3+  3+ 

and  YbFeO^.  The  Fe  <•  Yb  exchange  was  observed  in  the  optical  spectrum 

3+  “1 

of  Yb  in  YbFeO^.  The  exchange  splittings  were  large,  6  to  10  cm  , 

3+  3^ 

although  the  Yb  and  Fe  spins  are  at  nearly  right  angles,  indicating 
a  large  antisymmetric  exchange.  Further  studies  are  suggested  to 
illuminate  the  orthoferrite  magnetic  structure. 
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INVESTIGATION  OF  THE  ENERGY  LEVELS  OF  MAGNETIC  IONS 
IN  THE  COMPLEX  METAL  OXIDES.  ESPECIALLY  IN  THE  ORTHOFERRITES 

I.  INTRODUCTION. 

The  rare  earth  orthoferrites  are  a  family  of  mixed  oxides  of 

chemical  formula  mV^03  which  crystallize  in  a  slightly  distorted 

perovsklte^  ^  form.  M^  is  typically  a  trlvalent  rare  earth  ion  and  M^^ 

a  trlvalent  transition  metal  ion,  though  variations  from  this  typical 

9-13 

form  (involving,  say,  Ca  and  Ba  substitution)  have  been  studied. 

We  will,  in  this  report,  restrict  our  attention  primarily  to  the  cases 

where  M^  is  Y^^,  Gd^  or  Yb^,  and  M^^  is  Fe^,  Al^  or  Ga^. 

The  rare  earth  orthoferrites  are  interesting  magnetically 

14 

because  they  are  in  many  Instances  canted  antlferromagnets.  The 
magnetic  ions  group  into  sublattice  whose  magnetization  is  equal 
and  vdiose  direction  is  nearly,  but  not  quite  exactly,  antiparallel. 

There  is  therefore  a  small  permanent  magnetic  moment  due  to  the  failure 
of  the  sublattice  magnetizations  to  cancel  fully.  The  "canting  angle" 
is  the  angle  idiich  the  magnetic  sublattices  make  with  a  truly  antiparallel 
configuration,  and  the  origin,  temperature  dependence,  etc.  of  the  canting 
is  of  considerable  physical  interest. 

In  this  report  we  will  examine  the  spin  configurations  possible 
in  the  orthoferrites  from  both  a  macroscopic  and  microscopic  point  of 
view.  The  work  reported  is  the  initial  portion  of  the  work  aimed  at 
understanding  the  magnetic  properties  of  these  materials  from  an  atomic 
point  of  view. 
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II 


CRYSTAL  STRUCTURE  AMD  SPIN  CONFIGURATION  OF  ORTHOFERRITES. 


The  crystal  structure  of  the  rare  earth  orthoferrites  and  their 

related  compound,  the  orthoalumlnates,  orthogalllates,  orthochromites, 

etc.  have  been  studied  In  considerable  detail  by  Geller^  ^  and  associates 

0 

In  this  country  and  by  Bertaut  and  associates  In  France.  We  have  unit 
cell  dimensions  for  a  wide  range  of  combinations  but  have  at 

present  detailed  lattice  parameter  Information  -  e.g.  exact  oxygen 
positions  >  only  for  gadolinium  orthoferrite. ^ 

The  orthoferrites  and  their  Isomorphs  (we  will  use  the  term 
orthoferrite  In  this  report  as  the  generic  term  for  the  family) 
crystallize  In  the  cubic  perovsklte  structure  or  In  two  distorted 
modifications  of  the  perovsklte  structure,  one  orthorhombic  and  one 
rhombohedral.  The  compotonds  of  greatest  Interest  to  us,  the  Iron¬ 
bearing  rare  earth  orthoferrites,  all  crystallize  In  the  orthorhombic 
structure^  belonging  to  spac^  group  -  Dbnm.  There  are  four  distorted 
perovsklte  pseudo-cells  In  an  orthorhombic  unit  cell,  giving  four 
Inequlvalent  Iron  sites  and  four  Inequlvalent  rare  earth  sites  per  unit 
cell.  The  orthorhombic  unit  cell  Is  drawn  In  Fig.  1.  The  rare  earth  Ions 
have  been  omitted  In  Fig.  1  to  simplify  the  picture.  The  Iron  sites  and 
rare  earth  sites  of  the  unit  cell  are  shown  separately  In  Fig.  2  and  3, 
with  certain  additional  notations  which  will  be  utilized  later.  The  Iron 
positions  are  special  position  (1/2,  0,  0;  1/2,  0,  1/2;  0,  1/2,  0;  0,  1/2, 
1/2),  but  all  oxygens  and  rare  earth  Ions  are  displaced  by  as  much  as  a 
few  tenths  of  an  Angstrom  from  the  Idealized  sites  drawn  In  Figs.  1 
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through  3.  niese  distortions  will  be  discussed  later  relevant  to  the 
paramagnetic  resonance  and  optical  absorption  measurements. 

One  can  tell  a  great  deal  about  the  allowed  spin  configurations 
of  the  orthoferrite  on  the  basis  of  magnetic  symmetry  considerations.  The 
magnetic  groups  are  obtained  by  adding  to  the  symmetry  point  group 
elements  a  time  reversal  operation,  R,  which  reverses  all  angular  moments 
and  magnetic  fields. 

The  orthorhombic  space  group  to  which  the  orthoferrites 
belong  has  the  point  symmetry  group  elements 

Sx"  Sy'  ^2z'  **'y' 

where 

E  Is  the  Identity  operator 

^2y'  ^2z  rotations  of  n  about  the  designated  axes 

1  Is  the  Inversion  operator 

<s"^  are  reflection  syimnetry  operations. 

The  addition  of  R,  time  reversal,  to  the  point  syimnetry  elements  of  a 
crystal  containing  angular  momentum  results  In  the  generation  of  several 
magnetic  groups.  The  point  group  Itself  remains  a  symnetry  grotqi  and  a 
number  of  additional  groups  are  generated,  each  of  which  contains  R  in 
half  its  elements. 
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The  orthorhombic  group  D2j^  gives  rise,  in  the  presence  of  angular 
momenta,  to  the  following  magnetic  groups: 


(1) 

E 

I 

Sx 
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<S'z 
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^Z 

RI 

‘'"2y 

‘‘Sz 

R(5" 

(7) 
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^2y 

RI 

^‘=2z 

R(5" 

(8) 
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"2z 

<J'x 

<^y 

RI 

^‘'2x 

"‘=2y 

R<r 

We  now  look  at  the  orthoferrite  to  see  to  which  of  these  groups 
it  may  belong.  The  point  group  operations  transform  the  various  Ion  sites 
Into  one  another,  so  we  must  first  develop  a  scheme  of  labeling  of  the 
cation  sites  on  which  the  angular  momentum  (spin)  is  located.  We  will 
use  the  scheme  indicated  in  Figs.  2  and  3,  which  is  consistent  with  the 
notation  of  Heeger^^  in  his  discussion  of  the  Isostructural  KMnF^,  and  of 
Treves^^  in  his  discussion  of  magnetic  susceptibilities  in  the  orthoferrites. 
The  convention  is  slmlltu:  for  the  iron  and  for  the  rare  earth  sites: 
Horizontal  nearest  neighbors  (1,4)  (2,3) 

Vertical  nearest  neighbors  (1,2)  (3,4) 

Next-to-nearest  neighbors  (1,3)  (2,4). 

Table  ll-l  lists  the  possible  permutations  of  the  cation  sites  and  gives 
the  corresponding  symmetry  operations  which  produce  them 
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Permutation 

Fe 

R.E. 

Identity 

-  E,  I 

E, 

2-^3 

^2y'  ‘^k 

1  3,  2  <«->  4 

^2y’  'Sy 

^2y>  '^y 

la-^2,  3*-*>4 

Sz^  '^z 

Sz>  I 

TABLE  II- 1 


We  first  note  that  the  inversion  operator,  I,  carries  the  iron 
sites  back  into  themselves  (without  reversing  angular  momentum).  The 
operation  RI,  which  would  map  iron  sites  onto  themselves  and  then  reverse 
their  angular  momentum,  therefore  cannot  be  an  element  of  the  orthoferrite 
magnetic  symmetry  group.  We  therefore  eliminate  immediately  magnetic 
symmetry  groups  5  through  8  of  II-2  for  the  orthoferrite.  We  can  now 
deduce  the  spin  configurations  compatible  with  magnetic  symmetry  for  groups 
1  through  4,  since  we  know  both  what  the  group  operators  do  to  the 
individual  angular  momenta  and  also  what  permutations  are  Induced.  For 
Instance,  denoting  the  Fe  angular  momenta  by  S^(l  •>  1,  ...4.)  and  the 
rare  earth  angular  momentum  by  J^,  we  have 


C 

C 
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2x^lx 

2x®ly 

2x®lz 


"lx 

-S 


"4x 


ly  "4y 


'Iz 


’4z' 
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Therefore  any  group  containing  element  corresponds  to  a  spin 
configuration  with 

hx  "  hx’  hy  “  ■  %>  hz  ■  ■  hz  11-^ 


Similarly 


lx 

ly 

Iz 


'lx 


4x 


II-5 


so  spin  configurations  corresponding  to  a  magnetic  group  containing  (RC2  ) 
will  have 


^Ix  "  ■®4x  ^ly  "  hx  hz  "  hz- 

Following  this  procedure  one  obtains  the  magnetic  structure  corresponding 
to  each  group.  Before  listing  these  configurations  we  further  note  one 
operator  of  special  Interest,  (S~ which  transforms  the  rare  earth  sites 
into  themselves.  For  groups  containing 


^z  '^ix  -  -  •'ix  -  •'ix 

^  ■'ly  ■  ■  •’ly  ■  ■’l,  II-7 

'Q  ■'ll  ■  ■'i.  ■  •’lx- 


Hence 

•^ix  "  “^iy  *■  °  II-8 

in  groups  containing 
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Collecting  the  results  of  such  considerations  one  has  the  allowed 


configurations: 
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We  have  pursued  the  problem  thus  far  with  the  only  Information 
input  being  the  crystalline  symmetry  of  the  orthoferrite.  We  now  Inject  a 
physical  as suiq)t ion  concerning  the  magnetic  structure,  namely  that  the 
interaction  between  nearest  neighbor  iron  ions  be  antiferromagnetic, 

$2  and  are  then  constrained  to  be  essentially  antiparallel  to  and 
S^.  We  find  that  magnetic  group  4  does  not  allow  of  such  a  configuration 
so  we  drop  it  from  further  consideration  in  this  report  on  the  antiferro¬ 
magnetic  orthoferrites.  (It  is  however  the  appropriate  group  for  certain 

ferromagnetic  configurations  in  the  orthochromite  and  orthomanganate 
9  10 

family  '  ).  The  other  three  groups  give  one  pure  antiferromagnetic 

and  two  canted  antiferromagnetic  configurations  as  follows 


Magnet Ic 
Symnetry 

Group 

Predominant  Fe  spin 
direction 

Canting 

direction 

1 

y 

None 

2 

z 

X 

3 

X 

z 

TABLE  II- 2 

All  three  configurations  are  actually  found. Type  3  Is  the 
commonest  form;  it  is  the  spin  configuration  for  yttrium  orthoferrite  at 
all  temperatures,  and  of  nearly  all  the  rare  earth  orthoferrites  at  room 
tenq;>erature  and  above.  At  lower  temperatures,  usually  below  100°  K,  many 
of  the  rare  earth  orthoferrites  switch  to  configuration  (2)  and  a  few  to 
configuration  (1). 
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The  stable  configuration  Is  probably  determined  by  the  anisotropies 
of  the  separate  spin  systems.  If,  for  Instance,  the  local  crystal  field 
effects  were  such  that  the  Iron  Ions  In  YFeO^  occupied  their  lowest  energy 
position  when  their  spin  was  parallel  to  an  orthorhombic  a  direction 
(x  direction),  configuration  (3)  would  be  preferred.  If  a  y-dlrected  spin 
gave  lowest  energy,  the  antiferromagnetic  configuration  (1)  would  be 
preferred,  and  so  on.  Hie  paramagnetic  resonance  and  optical  absorption 
experiments  reported  below  were  Instituted  to  establish  these  correlations, 
and  the  relationship  of  these  spectra  to  the  actual  spin  configurations 
will  be  discussed  in  section  IV-C. 


III.  MAGNETIC  SUSCEPTIBILITIES  AND  ANTIFERROMAGNETIC  RESONANCE  MODES. 

Before  proceeding  to  the  atomic  spectroscopy  and  Its  Interpretation 
we  shall  first  examine  some  of  the  more  macroscopic  magnetic  properties 

of  the  orthoferrites.  It  Is  Indeed  from  such  macroscopic  magnetic  measure- 

17-21  9-13 

mentS)  plus  some  very  elegant  neutron  diffraction  experiments, 

that  we  have  most  of  our  present  Information  on  the  orthoferrites. 

An  extensive  set  of  magnetic  susceptibility  measurements  have 

19  20  17  21 

been  made  by  the  Grenoble  group  In  France  '  and  by  Bozorth  ’  and 
associates  at  Bell  Telephone  Laboratories, though  the  bulk  of  the  results 
have  never  been  published  in  detail.  The  principal  results  may  be 
summarized  as  follows. 

The  magnetization  of  the  rare  earth  ferrites  may  be  described 
by  a  permanent  moment  plus  a  susceptibility  (much  as  for  the  garnets) 

<r  »  <5^  +  XH.  III-l 

The  ferromagnetic  (or  permanent)  component  Is  about  .05  Bohr  magnetons 
per  molecule  at  room  temperature  for  all  the  rare  earth  orthoferrites.  The 
direction  of  <5"  ^  Is  always  in  either  an  orthorhombic  £  or  orthorhombic  a 
direction.  For  all  of  the  crystals  observed  except  SmFeO^  the  c  axis  Is 
the  direction  of  room  temperature.  For  an  orthoferrite  with  a 

magnetically  Inactive  rare  earth  (such  as  Y  or  La)  the  permanent  moment 
<S~ ^  remains  essentially  constant  In  magnitude  and  direction  down  to  liquid 
Helium  temperatures.  For  a  rare  earth  orthoferrite  containing  a  rare 
earth  such  as  Holmium,  Erbium,  or  Dysprosium  the  permanent  moment  Increases 
substantially  as  the  tenperature  Is  lowered,  and  usually  flips  to  the 
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a  direction  or  drops  suddenly  to  zero  at  a  transition  temperature.  Some 
exceptional  behaviors  appear;  gadolinium  with  Its  large  magnetic  moment 
appears  to  act  as  If  It  were  magnetically  Inert  In  the  orthoferrite. 

The  general  features  of  the  magnetic  behavior  we  can  understand 

on  the  basis  of  our  foregoing  magnetic  group  theory,  and  by  extrapolation 

22 

from  other  magnetic  materials,  notably  the  garnets. 

Let  us  first  assume  that  the  Iron  Ions  are  coupled  antlferro- 
magnetlcally  with  a  strong  Interaction  (T^  ^680°  K).  The  rare  earths 
are  then  coupled  to  the  Iron  Ions  with  a  rather  weak  exchange  (T^  fssi20°  K) . 
At  temperatures  appreciably  higher  than  V  such  as  room  temperature,  the 
rare  earths  are  negligibly  polarized  by  the  exchange  Interaction,  and  <S~^ 

Is  due  primarily  to  the  Iron  Ions.  We  can  account  for  the  observed 
properties  In  this  temperature  region  by  assuming  that  the  Iron  Ions 
occupy  the  spin  configuration  given  by  magnetic  group  3  of  Table  II-2. 

The  Iron  Ions  then  point  predominantly  In  the  plus  and  minus  a  (x) 
direction  but  are  canted  slightly  In  the  z  direction.  A  canting  of 
about  .01  radian  would  produce  the  moment  observed;  ,05  Bohr  magnetons 
per  molecule  out  of  a  possible  5.  As  the  temperature  Is  lowered  the 
weak  rare-earth  to  Iron  exchange  polarizes  the  rare  earths  more  and 
more  effectively,  giving  rise  to  an  Increasing  (S'^.  As  the  rare  earths 
become  more  polarized,  anisotropy  effects  arising  In  the  rare  earth  Ion 
lattice  become  more  effective,  and  the  whole  magnetic  structure  can  be 
flipped  to  that  of  magnetic  groups  No.  2  (Fe  spins  In  +  z  direction. 
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canted  toward  x)  or  to  the  antiferromagnetic  configuration  of  symmetry 
group  No.  1. 

The  general  explanation  raises  a  number  of  specific  questions. 

For  Instance^  why  do  the  Iron  sublattices  choose  to  cant  with  respect  to 

one  another?  An  Isotropic  exchange  would  cause  the  spins  to  minimize 

their  energy  by  being  rigorously  antiparallel.  The  canting  could  be 

caused  either  by  single  Ion  anisotropy,  or  by  antisymmetric  exchange. 

The  former  Is  caused  by  the  Ions  "feeling"  the  crystal  lattice  through 

crystalline  electric  field  effects;  the  second  Is  an  exchange  of  the  form 

D  •  X  $2  lAlch  tends  to  make  the  spins  be  perpendicular  to  each  other. 

23  24 

Dzlaloshlnskl  and  Morlya  have  shown  that  such  an  Interaction  Is  allowed 
(on  symmetry  grounds)  in  the  orthoferrite  and  have  suggested  an  atomic 
mechanism  for  the  interaction.  Treves^^  has  attempted  to  Identify  the 
canting  mechanism  on  the  basis  of  torque  measurements  and  concludes 
that  the  antisymmetric  exchange  Is  responsible.  We  have  come  to 
substantially  the  same  conclusion  based  on  other  arguments,  detailed 
subsequently. 

Let  us  first  derive  the  magnetic  susceptibility  and  magnetic 

resonant  modes  of  the  canted  antlferromagnet.  This  calculation  has 

been  reported  In  detail  In  a  scientific  report,  "Resonance  and  RF 

25 

Susceptibilities  In  Orthoferrites"  by  G.  F.  Herrmann,  and  we  will 
only  summarize  the  results  here. 

We  restrict  our  attention  first  to  the  simplest  case,  YFeO^, 
where  the  rare  earth  sites  are  occupied  by  non-magnetlc  ions.  The  appropriate 
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magnetic  group  la  No.  2,  the  irons  point  predominantly  In  the  x  direction 
but  are  canted  slightly  In  the  z  direction.  There  are  two  antiferromagnetic 


resonant  modes  of  the  system;  cu  ,  for  \dilch  the  precession  ellipse 

2 

of  each  sublattice  magnetization  Is  elongated  In  the  xz  plane;  and  co  , 

xy 

for  vdilch  precession  ellipse  Is  elongated  In  the  xy  plane. 

We  will  derive  the  required  properties  from  an  expression  for 
the  free  energy.  It  Is  convenient  to  express  the  energy  In  units  of 
magnetic  field.  Defining 

m-|«i|-|J!2|. 

idiere  and  ^  are  the  sublattice  magnetization  vectors  pointing  approxi¬ 
mately  In  the  plus  and  minus  x  directions  respectively^  we  shall  put 

V  -  W/M, 


where  W  Is  the  free  energy.  We  shall  then  write  W  as  a  function  of  the 
unit  vectors 

Rl  - 

R2  " 

R,  ■  M.  /M  and  R,  *  M„  /M, 
and  retain  only  quadratic  terms  in  R^^  and  R^. 

The  most  general  form  of  the  magnetic  anisotropy  tensor  for 
allowed  by  symmetry  considerations  Is  then 


■(A  X,  ^  +  A 
'  XX  1  zz 


2 

z,  +  A  X, 

1  xz  1  1 
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III-3 


III-4 


III-5 
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where  A _ .  A  ,  and  A  ere  macroecopic  energy  coefficients  related  to 

XX  zs  xc 

the  crystal  field  splittings.  Similarly  the  anisotropy  tensor  for  M2  Is 


■(A  ^  X, 
^  XX  2 


+  A 


zz 


-  V2>' 
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The  symnetrlc  and  antisymmetric  exchange  terms  are  given  by 

h 

and 

-D  (x^Z2  -  *2*1^ 
respectively. 

The  normalized  free  energy  of  the  system  Is  then 
V  -  E(Rj  •^)  -  D(XjZ2  -  X2Z^) 

-  *2^^  ”  ^zz  ^*1^  *2^^  III-IO 

"  <*1*1  -  V2)- 

One  then  solves  for  the  equilibrium'  position  of  and  H2>  their  normal 
modes  of  oscillation,  and  their  perturbation  by  an  external  magnetic  field. 

This  Is  done  In  the  report  referenced  above,  and  we  shall  list  the  Important 
conclusions  here. 

We  have  for  the  magnetic  susceptibilities 


(f  -  %  IL  +  2C  H  H 
”  X  XX  X  XXZ  X  z 

cT  "  %  H 

Oy  /'-yy  y 

(S’  ~  <r  +XH+C 

^  Z  ZO  '^ZZ  Z  XXZ  X 

where 
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Certain  approximations  have  been  made  In  the  above  expressions  which  amount 
substantially  to  the  assunqition  that  the  Isotropic  exchange  energy  Is  large 
compared  to  the  other  energies  involved^  and  that  the  canting  angles  are 
small. 

Since  a  measure  of  the  canting  effect  of  single  Ion 

anisotropy^  and  the  antisymmetric  exchange  parameter,  enter  the 

expressions  f or ^  ,  and  g-  to  first  order  In  the  same  form,  (A  +  D), 

XX  xz 

it  is  not  possible  to  determine  these  parameters  separately  from  susceptibility 

measurements  until  one  Invokes  a  measurement  of  C  ,  a  very  small  higher 

xxz 

order  coefficient.  Susceptibility  measurements  are  therefore  a  poor  tool 
for  distinguishing  these  two  effects  analytically.  However,  the 

zo 

measured  Implies  that 
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or  that  the  "canting  field"  be  approximately  one  one-hundreth  of  the  exchange 
field.  Since  the  equivalent  exchange  field  Is  on  the  order  of  10^  oe.,  we 


15 


deduce  that  the  "canting  field,"  a  fictitious  field  at  right  angles  to 

the  exchange  field  capable  of  producing  the  observed  canting,  must  be  on 

the  order  of  10^  oe.,  which  is  one  or  two  orders  of  magnitude  greater  than 

3+ 

the  anisotropy  fields  usually  associated  with  the  Fe  Ions.  (A  more 
quantitative  discussion  of  this  Is  given  below).  We  are  therefore  Inclined 
by  the  susceptibility  measurements  to  suspect  antisymmetric  exchange  as 
the  source  of  canting. 

A  better  technique  for  separating  out  single  Ion  anisotropy 

effects  from  antisymmetric  exchange  effects  Is  offered  by  antiferromagnetic 

resonance.  The  two  sublattice  system  has  two  resonance  modes  of 

frequency  <u  and  <u  as  mentioned  above  and  described  In  the  referenced 
z 

report.  In  the  absence  of  an  external  field  the  resonant  frequencies 
are 


co^^0)/y^  ■ 

co^,^0)/y2 


4E(A  -  A  )+4A  (A  -  A  )+4A  (A  +  D) 

XX  ZZ  XX  XX  zz  xz  xz 

■  4E(A  )+4A  (A  -  a  )+(A  +  D)^. 

XX  XX  XX  zz  xz 
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The  dominant  term  for  both  modes  is  the  typical  /H  H  type  term.  For  the 
7  3  5 

orthoferrite  fisilO  and  ?%4l0  — e-lO  (guesses  based  on  T^  and  typical 
3+ 

Fe  anisotropy  fields)  so  the  resonant  frequency  might  fall  at  frequencies 

ranging  from  y  x  10^  (or  280  kMc)  to  y  ^  3  x  10^  (or  1,000  kMc).  On  the 

other  hand,  the  (A  -  A  )  cancellation  In  the  co  mode  might  bring  the 
’  '  XX  zz  z 

frequency  downward  to  100  kMc  or  less.  Magnetic  resonances  of  some  sort 

26 

were  once  seen  In  the  orthoferrites  by  Dr.  Simon  Foner,  at  the  Lincoln 
Laboratories,  at  70  kMc  In  magnetic  fields  ranging  from  a  few  kllogauss  up 
to  120  kllogauss. 
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In  the  presence  of  magnetic  fields  the  resonance  frequencies 

become; 

In  z  direction  (£  axis) 

-  a3^^(0)/Y^  +  (A^^+  D)H^  111-15 

cu  -  0)  ^(0)/y^  +  (5A  +  D)H  +  H 

xy  '  xz  '  xz  o  o 

Expressions  for  with  In  other  directions  are  substantially  more 

conqillcated  and  are  given  In  the  referenced  report.  One  sees,  however, 
that  from  the  measurement  of  oo  and  cu  as  a  function  of  H  ,  the  single 
Ion  canting  effect  (measured  by  A  )  and  the  antisymmetric  exchange  effect 

XZ 

(measured  by  D)  could  be  analytically  separated. 

We  conducted  a  search  for  the  antiferromagnetic  resonance  modes 
at  frequencies  from  55  to  65  kMc  and  In  magnetic  field  up  to  15  kllogauss. 

No  resonances  were  observed.  We  conclude  that  our  observation  frequencies 
probably  were  not  high  enough.  The  closing  of  our  laboratory  precludes 
the  extension  of  this  search  to  higher  frequencies  within  the  present 
contract. 
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IV.  PARAMAGNETIC  RESONANCE:  IRON 

A.  Experimental. 

Yttrium  orthoferrite  Is  the  simplest  interesting  orthoferrite, 

In  the  sense  that  It  has  only  one  kind  of  magnetically  active  Ion  present, 

trlvalent  Iron.  The  magnetic  properties  of  the  yttrium  orthoferrite  are 

3+ 

presumably  fully  determined  by  the  nature  of  the  ground  state  of  the  Fe 

Ion,  by  Its  Interaction  with  Its  crystalline  environment,  and  by  its 

3+ 

Inter action  with  other  Fe  Ions. 

3+ 

The  properties  of  the  Individual  Fe  Ion  and  Its  interaction  with 

Its  crystalline  environment  can  be  determined  through  paramagnetic  resonance 

3+  6 

spectroscopy.  The  ground  state  of  the  Fe  Ion  is  a  spectroscopic  S  state, 

with  L  0  and  S  »  5/2.  Since  It  Is  an  S  state  Its  charge  cloud  Is 

essentially  spherically  symmetric,  and  Its  energy  of  Interaction  with  the 

crystalline  electric  potentials  small,  typically  a  few  tenths  of  a  wave 

number,  (or  a  few  tens  of  kilomegacycles).  To  study  the  isolated  Ion  It 

must  be  substituted  dllutely  Into  a  diamagnetic  Isomorph  of  the  crystal 

of  magnetic  Interest.  We  have  used  YAIO^,  yttrium  orthoaluminate,  for  our 
27 

host  crystal  In  the  paramagnetic  and  optical  studies.  YGaO^  would  have 
been  a  more  Ideal  diamagnetic  host  since  trlvalent  gallium  and  trlvalent 
Iron  have  nearly  the  same  radius,  whereas  aluminum  Is  substantially  smaller. 
However,  YGaO^  has  never  been  reported  as  a  stable  phase,  and  all  our 
attempts  to  grow  YGaO^  were  unsuccessful  (resulting  usually  In  the 
garnet  structure 
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A  substantial  portion  of  the  report  period  was  In  fact  devoted 
to  the  problem  of  growing  single  untwlnned  crystals  of  YAIO^.  After  some 
considerable  exploratory  work  It  was  found  that  a  melt  containing  2  mole 
percent  Al203>  2  mole  percent  Y202>  88  mole  percent  PbO  and  8  mole  percent 
PbF2  could  be  made  to  yield  regularly  crystals  approximately  2  nan  on  a 
side.  However^  these  crystals  were  found  to  be  almost  invariably  multiply 
twinned.  Since  small  additions  of  SIO2  had  been  previously  found  to 
affect  substantially  the  growth  habit  of  the  aluminum  garnets,  a  small 
addition  of  SIO2  to  the  melt  was  tried,  and  an  addition  of  approximately 
1  mole  percent  of  SIO2  was  found  to  suppress  the  twinning  strongly. 
Spectrographlc  analysis  has  shown  that  the  SI  does  not  enter  the  ortho* 
alumlnate  crystal  In  the  process  of  modifying  the  growth  habit;  the  SI 
content  of  YAIO^  grown  from  SIO2  doped  melts  and  undoped  melts  was  sub¬ 
stantially  the  same.  The  crystals  used  throughout  the  paramagnetic 

resonance  studies  described  below  were  grown  from  the  melt  composition 

3+  3^  3+ 

described  above,  plus  appropriate  Fe  ,  Gd  ,  or  Yb  dopings. 

The  paramagnetic  resonance  studies  were  performed  primarily  on 

a  35  kMc  superheterodyne,  magnetic  resonance  spectrometer.  Most  of  the 

data  were  taken  at  liquid  nitrogen  temperatures,  with  a  few  measurements 
o  3+ 

at  4.2  K.  The  Fe  site  In  the  orthoferrite  Is  of  low  symmetry,  having 
only  the  Inversion  symmetry.  Since  there  are  four  magnetically  Inequlvalent 
Iron  sites  per  unit  cell,  and  five  transitions  (of  ■■  1  type)  per 
Ion,  one  expects  the  principal  spectrum  to  consist  of  20  absorptions 
for  a  general  orientation  of  the  external  magnetic  field,  H^.  For 
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In  the  principal  planes  of  the  crystal  (e.g.  the  ab  plane)  the  number  of 

magnetically  Inequlvalent  sites  Is  reduced  to  two,  and  along  the  b,  or  c 

axes  the  Iron  sites  are  all  magnetically  equivalent. 

The  Iron  In  the  orthoalumlnate  Is  octahedrally  coordinated,  that 

Is, has  six  oxygen  nearest  neighbors  In  a  nearly  cubical  array.  The  oxygens 

are  all  displaced  from  their  true  cubic  positions  by  amounts  of  a  few 
®  6 

tenths  of  an  Angstrom  out  of  an  Ion  separation  of  approximately  two 
o 

Angstroms.  The  position  of  Ions  Is  described  more  fully  In  the  crystal 
field  calculation  of  section  IV-C.  For  present  purposes  It  Is  sufficient  to 
note  that  the  low  symmetry  of  the  site  means  that  one  has  no  a  priori 
knowledge  of  the  orientation  of  the  principal  axes  ot  the  crystal  field 
distortions,  and  indeed  no  guarantee  that  there  would  be  any  slnq>le 
relationship  between  the  principal  axes  of  the  quadratic  potential  and 
those  of  the  biquadratic  (fourth  power)  potential. 

3+ 

The  paramagnetic  resonance  absorption  spectrum  of  the  Fe  Ion 

in  YAIO-  was  first  recorded  for  H  in  the  a-b  plane  and  for  H  in  the  b-c 
j  o  ■  o 

plane.  These  data  are  presented  In  Figs.  4  and  5.  It  is  not  obvious 
from  these  data  that  any  principal  axes  were  intercepted.  However,  the 
extremal  spaclngs  of  the  outermost  transitions  at  about  25°  from  the  b 
axis  In  the  a-b  plane,  and  the  near  symmetry  In  the  spaclngs  of  a  group 
of  five  lines  at  this  orientation  (the  five  transitions  Identified  as 
sites  1  and  2  In  Fig.  4),  suggested  that  this  direction  was  nearly  a 
principal  direction  of  the  dominant  perturbation.  We  expected  the  dominant 
perturbation  to  be  quadratic  In  the  spin  conponents,  as  Is  usually  the 
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3+ 

case  for  the  Fe  spectrum.  We  therefore  took  data  with  In  a  plane 
containing  the  £  axis  and  the  extremum  In  the  a-b  plane.  These  data 
showed  that  the  extremum  was  located  In,  or  within  two  or  three  degrees 
of,  the  a-b  plane.  Concluding  then  that  this  direction,  25°  from  b  In 
the  a-b  plane,  was  a  principal  axis  of  the  quadratic  energy  tensor,  we 
took  data  with  In  the  plane  perpendicular  to  this  direction  (hence¬ 
forth  called  direction  3)  and  obtained  the  spectrum  shown  In  Fig.  6. 

B.  Determination  of  the  Spin  Hamiltonian. 

The  fitting  of  the  spin  Hamiltonian  for  a  spin  5/2  or  7/2 

particle  to  the  experimental  data  for  a  low  synnetry  site  presents  a 

formidable  problem.  In  the  present  Instance  the  Iron  sice  has  only 

Inversion  synmetry;  we  can  from  this  Information  eliminate  only  terms  in 

the  potential  which  correspond  to  Y^(0^)  where  m  Is  odd.  We  have 
0  4-2  0  4-2  4-  4 

therefore  Yj  ,  Y2-  ,  Y^  ,  Y^-  and  Y^-  terms  remaining,  and  no  5  priori 
knowledge  of  the  principal  axes  of  the  quadratic  terms  or  the  relationship 
these  axes  bear  to  the  axes  In  \dilch  the  fourth  order  potential  takes  Its 
simplest  form.  Each  guess  at  the  eight  constants  In  the  spin  Hamiltonian 
leads  to  a  6  X  6  matrix  for  each  site,  to  be  diagonalized  at  each  H^ 
orientation  for  comparison  with  the  experimental  data.  The  labor  Involved 
In  such  a  trial  and  error  fitting  Is  enormous,  even  with  Inspired  guessing, 
and  convergence  of  the  procedure  Is  In  no  way  guaranteed. 

It  was  necessary,  therefore,  to  develop  an  analytic  attack  which 
circumvented  the  necessity  for  repeated  dlagonallzatlon  of  the  energy 
matrix.  This  attack  proceeds  as  follows  and  on  the  principles  outlined 
In  the  next  paragraphs. 


The  grotind  state  energy  level  of  a  paramagnetic  S  state  Ion 


in  an  external  field  Is  customarily  given  In  terms  of  a  spin  Hamiltonian 
of  the  form 

gp  H  •  S'+ 

The  first  term  gives  the  Zeeman  energy  of  Interaction  with  the  external 
field,  and  ^^(S)  Is  the  anisotropy  energy  which  depends  essentially  on 
the  orientation  of  the  spin  vector  ?  with  respect  to  the  crystalline  axes. 
For  an  S  state  Ion,  g  can  be  assumed  to  be  isotropic  and  have  the  value  2. 
Hie  procedure  below  could  be  generalized  to  an  anisotropic  g  but  there  Is 
no  merit  for  our  present  purposes  In  so  doing. 

One  wishes  to  determine  ^  (S)  from  the  paramagnetic  resonance 

& 

experiments.  If  the  Zeeman  energy  of  g^H  •  S  is  sufficiently  large, 
all  off-diagonal  elements  of  ^  will  have  negligible  effect  on  the 
spectrum,  and  can  be  deduced  from  consideration  of  the  diagonal 
matrix  elements  only,  without  further  attention  to  the  off  diagonal 
terms.  In  the  present  case,  however, ^  ^  Is  not  negligible  with  respect 
to  g^H  so  a  more  elaborate  analytic  technique  is  Involved.  We  shall  show, 
however,  that  It  Is  possible,  from  symmetry  considerations,  to  process  out 
of  the  data  the  effect  of  the  off-diagonal  elements,  and  again  extract 
the  full  properties  of  from  the  diagonal  elements  alone. 

It  is  convenient  to  define  from  the  outset  two  sets  of 
coordinate  systems;  the  first,  denoted  by  subscripts  1-,  2  and  3,  being 
fixed  In  the  crystal  with  respect  to  the  given  Ion  site,  and  the  second. 


IV- 1 
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denoted  by  x,  y,  z,  such  that  z  coincides  with  the  applied  field  H. 

Is  usually  given  as  a  power  series  In  the  spin  components  $2  and 
S^.  We  shall  put 

e 

where  are  polynomials  of  order  ^  In  the  spin  components,  and  the 
factor  has  been  Included  to  give  In  units  of  magnetic  field,  as 
Is  connon« 


For  the  purpose  of  our  analysis  it  Is  Important  to  construct 

28“33  /) 

the  as  operator  equivalents  of  harmon'^c  polynomials  of  order 

In  the  directional  cosines  of  some  suitable  radius  vector.  This  can 
always  be  accomplished  by  the  addition  of  certain  spherically  symmetric 
terms  to  any  of  the  polynomials  or  their  operator  equivalents.  Once 
written  In  this  form,  can  be  expressed  as  a  sum  of  spherical  harmonic 
functions  Yj  In  any  coordinate  system,  since  for  any  fixed the  ^ 
transform  only  among  themselves  In  any  orthogonal  coordinate  trans¬ 
formation.  In  the  paramagnetic  case  the  Vj^  are  real  and  symmetric  In 
the  variables.  For  any  given  value  of  S  one  need  consider  only  terms 
with  ^  <  2S. 


We  now  apply  2nd  order  perturbation  theory,  assuming,^  small 

CL 

compared  to  g^H.  Quantizing  along  the  z  direction  one  can  write  the 
following  expression  for  the  energy  levels: 

2 


>  + 


M. 


(Mg-  Mg')gpH 


IV- 2 


IV- 3 
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Consider  first  the  diagonal  term  <  >.  We  expand  each  polynomial 

in  terms  of  operator  equivalents  of  the  spherical  harmonics  in  the 
(x,y,z)  system 

m  y  2 

Since  this  expansion  involves  a  transformation  from  the  (123)  to  the  (xyz) 
system^  the  coefficients  are  functions  of  the  orientation  of  H  or 

-  aJ  (0^,02,03)  ] 

where  ^represents  a  unit  vector  along  7.  The  diagonal  matrix  elements 
of  are  given  by 


“a  >  “  ^  ^  (S)  I  Mg  > 


IV- 6 


since  for  m  0  (S)  has  no  diagonal  elements.  To  obtain  the  angular 

dependence  of  ,  one  substitutes  the  directional  cosines  of  iT  for 

in  Eq.  (IV-5)  .  In  the  (xyz)  system  one  has,  by  definition,  a  -  a  »  0 

X  y  ' 

*  1.  Therefore,  since,  (0,0,1)  »  0  for  all  m  0 


y|  (0,0,1) 


IV- 7 


*  When  substituting  for  and  vice  versa,  one  must  replace  any 

polynomial  in  by  its  proper  operator  equivalent.  In  this  sense,  any 
equation  satisfied  by  is  also  satisfied  by  and  vice  versa. 
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and  (IV-6)  becomes 


(0^0^03)  <  M  J  I  > 


(0,0,1) 


The  angular  dependence  of  the  diagonal  matrix  elements  of  Is  thus  given, 

apart  from  numerical  factors,  by  the  "classical"  expression  (aj^a2a3). 

The  spectrum  Is  thus  the  product  of  two  parts.  One,  <  H  /  Y?  /  M  >, 

s  JL  s 

gives  the  fixed  structure  of  the  spectrum  as  determined  by  the  quantum 
numbers.  The  other  Is  a  scaling  factor  (a^a2Q!3)  which  depends  only  on 
the  orientation  of  H. 

We  give  In  Table  IV-1  a  list  of  a  few  of  the  more  useful 
polynomials  and  their  operator  equivalents;  consider  for  example  a  spin 
Hamiltonian  term  of  the  form 


V^(s)  =  D  [83^  -  3  8(8+1)]  +  E(8^^-  82^) 


one  then  has 


=  0(03^  -  3)  +  £(0^“-  OL^) 


IV- 10 


and,  making  use  of  Table  lV-1  and  Eq.  (IV-8), 

<  ^  I  ^"2  I  “s  >  -  2  -  3>  +  “3^>]  * 

[m^2  .  i  s(8+l)]  . 

The  problem  of  flndiig  the  spin  Hamiltonian  has  thus  been 
reduced  to  finding  the  orientation  dependence  of  <  |  >. 


IV- 11 
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o  2  1 

-  “3  -  I 

<  “s  I  Vi  “s  >  "  ■  3 

"  ih  V  ■  ^°“3^  +  ^ 

<  Mg  /  Y^°|  Mg  >  -  [35  Mg^  -  30S  (S+l)Mg^+  25Mg^-6S(S+l)+3S^(S+l) 

“  2^0  ^3^  -  315  03^  +  105  03^-  5] 

<  “s  i  VI  >  "  2^0  [2^^  ^8^  ■  315S(S+l)MgV  735Mg^  + 

105S^(S+l)^Mg^-  525S(S+l)Mg^+294Mg^-5S^(S+l)^ 

+  40S^(S+1)^-  6S(&fl)J 

''cub<^>  "  i  +  “2''  +  V  ‘  5  ] 

“  6  [h^  *  V  V  “  5  3S-1)]  . 

Also 

Y2° (0,0,1)  .  I 

Y^°(0,0,1)  =  ^ 

V(0,0,l) 


TABLE  IV- 1 

© 
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The  next  step  involves  the  extraction  of  <  M  |  |  M  >  from  the 

experimental  data,  by  eliminating  the  contributions  due  to  off-diagonal 
matrix  elements.  This  is  made  possible  by  the  following  symmetry 
considerations.  Suppose  that  we  reverse  the  signs  of  and  M^'  in 
Eq.  (IV-3).  Since  <  M  I -5^  I  M  >  and  |  <  M  1^  I  M  '  >  |  ^  do  not  change 
sign,  one  finds 


W(-M  -  gPHM  +<M^  /X  /M  >  - 


M  * 
8 


(M^  -  M^')gPH 


and,  hence,  the  sum  W(M  )-fW(-M  )  is  in  first  order  Independent  of  the  off- 

8  8 

diagonal  terms.  By  using  proper  linear  combinations  of  the  transition 
frequencies  one  can  thus  eliminate  the  off-diagonal  contribution  and 
determine  separately  the  orientation  dependence  of  each  V £  ,  For  any 
value  of  S,  there  are  always  exactly  enough  transitions  to  make  such  a 
determination  coiq>lete. 

3+ 

Let  us  now  turn  explicitly  to  the  solution  of  the  Fe  spectrum. 


For  Fe  ;  ^  ■  2  and  y( ^  is  given  by 


IV- 12 


=  gPV2(T)  +  g3V^(S). 


IV- 13 


The  energy  levels  can  be  written  out  according  to  Eq.  IV-3.  The 
diagonal  part  <  j X >  is  found  from  Eq.  IV-8,  using  values  for 
<  M  1  Y_°|  M  >  and  <  M  \Y, M  >  which  can  be  found  in  the  literature. 
According  to  standard  practice,  each  transition  is  characterized  in  terms 
of  the  magnetic  field  for  which  it  is  observed  at  a  given  fixed  frequency  v^. 
The  five  transitions.  Mg  “  '  2  ’  2^  '  2~^  '  2>  ‘  2~^h  2~*h 
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occur  at  magnetic  field  values  denoted  respectively  by  H2;  and 

One  obtains 

”1  “  “o  +  ^  1 

Hj  =  +  2k  °  -  I  A^°  +  ^2 

H-  =  H  +A,  IV-14 

i  o  i 

-  JAj"’  +  I  +4^ 

"5  ■  "o  -  **2°  •  V+^5 

where  «  hv^/g3>  and  where  represent  contributions  arising  from  the 
off-dlagonal  matrix  elements  of  ^  ,  From  IV-3  and  IV-12  It  follows  that 
*  4^  and  *  A^,  provided  Is  substituted  for  H  In  the 
denominator  of  the  off-dlagonal  contribution  In  the  above  equations. 

(This  substitution  In  fact  tends  to  Improve  the  approximation.)  One  can 
now  solve  (IV- 14)  to  give 

^2°  =  -h  &<»1-  “5>  +  ^  <»2-  »4>] 

IV- 15 

V  ■  K  -  "5  -  ^<“2  -  »4>]  • 

3+ 

We  now  apply  this  analysis  specifically  to  the  spectrum  of  Fe 
In  the  orthoalumlnate  as  plotted  In  Figs.  4,  5  and  6.  First  one  decides 
on  the  basis  of  Intuition  and  general  morphology  which  five  transitions  or 
groups  of  five  transitions  belong  to  the  same  site  for  any  one  orientation 
of  The  grouping  Is  Indicated  on  the  Figs.  4,  5  and  6.  One  then  extracts 
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the  parameter  and  for  each  orientation  of  the  magnetic  field.  These 
are  plotted  In  Figs.  7  and  8.  There  Is  one  set  of  axes  1^  2,  3  and  two 
crystal  field  parameters,  B  and  E,  such  that  all  the  points  on  the  curve 
can  be  fit  to  the  expression 


-  3/2  [0(03^-  3)  +  E(a^^-  02^)]  . 


IV- 16 


A  best  fit  to  the  data  was  obtained  with  D  •>  -  1.45  kilo-oersted,  E  ■>  0.30 


kilo-oersted,  and  the  1,2,3  axes  (for  a  particular  site)  are  defined  by 
the  following  transformation  upon  the  orthorhombic  axes,  a, b,c 


/ 


L 

c 

V  sin  22 


22°  Bin 

38° 

sin  22° 

22°  COB 

38° 

-sin  22° 

0 

CM 

CM 

cos  22° 

-.570 

.231 

.788 

.730 

-.295 

.616 

.375 

.927 

0 

3° 

r.O 


COS  38 
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The  angular  variation  of  A^^  In  each  of  the  experimental  configurations  was 
calculated  using  the  above  parameters  and  axes,  and  Is  plotted  as  the 
theoretical  curve  In  Fig.  7.  The  closeness  of  the  fit  Indicates  the 
accuracy  of  the  best  fit  solution. 

The  term  In  the  Hamiltonian  quadratic  In  the  spins  is  therefore 


-  -1.45  -  ^S(S+1)]  +  0.30  -  82^  J  IV- 18 

*  For  the  purpose  of  this  transformation  we  regard  1,  2,  3  and  a,b,c  as  unit 
vectors  In  their  respective  directions. 
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where  Sj^,  S2  and  are  the  spin  components  In  directions  1,  2  and  3  above. 
Expressed  somewhat  differently  for  subsequent  coiiq>ari8on  with  the  crystal 
field  calculation^  V2(^  may  be  written 

V2(S)  -  +  0282^  +  IV-19 

with  Dj^  “  +1.18,  D2  =  0.27  and  “  “1.45  kilo-oersted. 

We  turn  now  to  the  term  biquadratic  in  the  spins,  V^. 

This  term,  in  its  full  generality,  contains  nine  independent  parameters, 
and  a  complete  determination  of  all  these  parameters  from  the  spectrtim 
is  not  feasible.  We  will  make  some  8lnq>llfying  assumptions  concerning 
the  fourth  order  terms  and  let  the  fit  resulting  be  a  test  of  the 
validity  of  the  assunptions. 

The  simplifying  assumption  consists  of  assuming  that  has 
cubic  symmetry.  The  Fe^"*"  is  indeed  surrounded  by  6  0^"  ions  which 
form  a  slightly  distorted  octahedron  as  is  shown  in  Fig.  9.  The 
coordination  octahedron  is  tilted  and  rotated  by  substantial  angles 
(circa  15  )  from  the  coordination  octahedron  of  the  cubic  perovsklte, 
but  the  radius  vectors  to  the  0  are  still  nearly  orthogonal  to  one 
another  and  of  nearly  equal  length  (c.f.  equations  VI-1,  section  VI-A) 

We  therefore  take  V^(T)  to  be  of  the  form 

V4(S)  -  -  I  [s^  +  ^  S(S+1) (38^+38-1)]  IV-20 

in  the  appropriate  coordinate  frame  f  >  ^  •  It  then  follows,  from 

IV- 7  that 
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A^°  (a^,  a^,  a^)  -  [a^ f  “V  '  5 1 "  a(l-5^) 


where 


a^ 


0!*,  +  oty,  a-f  +  0!^  a 


IV-21 


IV-22 


Next,  one  attempts  to  fit  the  data  points  of  Fig.  8  with  a  value  of  a  and 

a  choice  of  axes.  The  theoretical  curves  of  Fig.  8  were  plotted  for 

a  -  -.035  Koe.  and  a  coordinate  frame  ^  ^  obtained  from  the  ortho- 

* 

rhombic  axes  by  the  transformation 


.807 

-.574 

.122 

.565 

.819 

.122 

-.174 

-.052 

.985 

IV-23 


The  generally  good  fit  can  be  taken  as  a  confirmation  of  the  approximation 
and  of  the  fit  parameters. 

It  Is  interesting  to  conq>are  the  axes  of  the  cubic  potential 
with  the  radius  vectors,  r^,  to  the  ligating  oxygens.  To  do  this  we 


IV-23  with 

* 

the  direction  cosines  of 

the  r^: 

/  .772 

-.620 

.187 

-  /  .577 

.794 

.196 

\  ^1  / 

^  -.139 

-.088 

.990 

IV-24 


One  sees  that  the  axes  of  the  cubic  potential  do  Indeed  point  essentially 
toward  the  oxygen  neighbors. 


*  Ke  again  regard  ^  ^  ,  'f  ,  r^,  and  a,b,c  as  unit  vectors  here. 
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The  sign  of  a  deduced  here  is  negative;  its  sign  relative  to  D 

can  be  deduced  from  the  spectrum  and  the  sign  of  D  can  be  deduced  from 

the  temperature  dependence  of  absorption  intensities.  The  negative 

34 

sign  of  a  is  in  contradiction  with  Geschwlnd's  general  conclusion 

3+ 

on  the  sign  of  a  for  Fe  in  octahedral  coordination,  but  the  experi¬ 
mental  evidence  seems  conclusive  in  this  case. 
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V.  PARAHAGNETIC  RESONANCE:  GADOLINIUM. 

A.  Experimental. 

3+ 

The  paramagnetic  resonance  experiments  on  Gd  In  YAIO^  were 

3+ 

performed  In  substantially  the  same  manner  as  those  on  Fe  in  YAIO^. 

The  spectrometer  frequency  was  approximately  35  kMc,  and  most  of  the 

experiments  were  done  at  liquid  nitrogen  tenq>erature8. 

The  rare  earth  ion  site  in  the  orthoferrite  is  in  a  mirror 

plane  perpendicular  to  the  c  axis.  Therefore  one  principal  axis  of  the 

relevant  potential  will  always  be  the  c  axls^  and  the  other  axes  are 

confined  to  the  a-b  plane.  Further^  the  rare  earth  sites  are  coupled 

in  pairs  by  an  inversion  through  the  iron  sites,  so  there  are  only 

two  magnetically  inequivalent  rare  earth  sites  for  the  general  orientation 

of  H  .  For  H  in  the  b-c  or  a-c  planes  the  four  sites  are  all 
0  o 

3+ 

magnetically  equivalent.  The  Gd  spectrum  is  therefore  relatively  more 
3+ 

simple  than  the  Fe  spectrum,  even  though  there  is  a  larger  number  of 

transitions  per  site  (seven  transitions  of  ±  1)  because  of  the 

3+ 

larger  spin,  S  >  7/2,  of  the  Gd  ion. 

3+ 

The  paramagnetic  resonance  absorption  spectrum  of  Gd  with 
in  the  a-b  plane  is  plotted  in  Fig.  9.  The  absorption  spectrum 
with  in  the  b-c  plane  is  plotted  in  Fig.  10.  All  necessary  information 
is  contained  in  these  two  spectra. 

B.  Determination  of  the  Spin  Hamiltonian. 

The  same  general  procedures  and  philosophy  that  were  used  in 
3+ 

solving  the  Fe  spectrum  will  be  followed  again  in  the  resolving  of  the 
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Gd  spectrum.  The  angular  behavior  of  the  diagonal  elements  of  the  spin 
Hamiltonian  will  be  used  to  derive  Information  on  the  full  spin  Hamiltonian, 
Including  off-dlagonal  elements. 

Since  the  spin  of  Gd  In  7/2,  one  must  Include  also  a  term  Vg(S} 

In  the  spin  Hamiltonian.  There  will  now  be  seven  transitions 

M  .5_..3  .3_^.l  .1_^1  1_^3  3^5  5_^7 

s  '  2’  2  2^  2  2^  2  2^  2  2^  2  2'  2  2’ 

which  occur  respectively  at  the  magnetic  field  values  H^,  H2  ,  H3,  H^,  H5, 

Hg,  H^.  One  finds: 

Hi  -  Ho  +  +  10  +  3  Ag°  +  ^ 

H2  -  H^  +  l*k°  +  5  A^°  -  7  A^°  +  A  2 

H,  -  H  +  2k  °  “  6A  °  +  7A,°  +  A,  V-1 

J  o  z  4  o  J 

-  H^  +  A^ 

H5  -  Hq  -  2k^  +  bk°  +  lk°  +  A  3 

«6  “  «o  ■  +  V  +  +  Ag 

-  Hq  -  6A2°  -  10A^°  -  3Ag°  +  A7  . 

Again,  for  the  off-dlagonal  contributions  one  has,  to  first  order, 

^  A  2  ~  A3  ~  solve  the  equations  to  give 

V  -  Tea  &<«r  «7>+®  (»2-  »6>  +  ^^«3“ 

V  “  3M  “7>'®  ^‘*2’  “6>  “  **5>]  ^’2 

V  “  ife  C<«r  «7>-'^  <«2-  +  5(H3-  H5)]  . 
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The  spectra  of  Figs.  9  and  10  were  processed  to  yield  the 
angular  forms  of  and  Ag°  In  the  and  b^c  planes  and  the  reduced 

data  for  k^  and  are  plotted  In  Figs.  11  and  12  respectively.  The 
variation  and  magnitudes  of  the  k^  values  were  within  the  experimental 
scatter  of  the  data,  so  the  k^  Information  was  not  regarded  as 
meaningful . 

We  have  taken  the  c  axis  to  be  the  z  or  distinct  axis  for  the 
various  terms  In  the  spin  Hamiltonian,  since  the  synmetry  of  the  site 
makes  this  a  unique  axis.  We  can  then  fit  the  k^  angular  dependence  to 


where  D  ■  +0.10  Koe. 

E  -  -0.35  Koe. 

The  principal  axes  1  and  2  are  derived  from  the  crystallographic  a  and 
b  axes  by  a  clockwise  rotation  of  28^  about  the  -tc  direction.  This  set 
of  principal  axes  bears  no  obvious  relationship  to  the  local  distortions 
at  the  rare  earth  site,  an  unusual  site  of  twelve-fold  oxygen  coordination. 

The  second  order  terms  may  be  also  written  In  less  conventional 
but  perhaps  more  Informative  form  as 


where  the  D. 


+  D2S2 


have  the  values 

-  -0.57  Koe. 

-  +0.47  Koe. 

-  +0.10  Koe. 


+  D  S 
c  c 


V-5 


V-6 
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If  we  turn  our  attention  now  to  the  biquadratic  term  In  the  spin 


Hamiltonian  we  find  a  much  more  conpllcated  situation  than  we  found 
34* 

for  the  Fe  spectrum. 

The  planes  In  which  data  were  taken  are  principal  planes  of 
the  coordination  cube  of  the  undlstorted  site.  If  the  site  syinnetry  remained 
essentially  cubic  In  the  actual  crystal,  the  variation  of  with  9 
could  be  represented  entirely  by  a  Cos (49)  function  In  both  planes. 

Inspection  of  the  data  on  A^°(9)  shows  this  to  be  far  from  the  actual 
case;  the  dominant  9  dependence  In  both  planes  Is  a  cos(29)  dependence. 
Therefore  we  cannot  make  the  simplifying  assumption  that  the  fourth  order 
potential  has  cubic  symmetry  and  that  one  field  magnitude  parameter  Is 
sufficient  to  characterize  the  potential.  Fortunately,  the  nine  general 
coefficients  of  the  potential  are  reduced  to  five  by  the  mirror  plane 
symmetry  of  the  rare  earth  site  (all  are  eliminated  for  idilch  m 

Is  odd).  can  therefore  be  written 

V^(^)  -  A(35a3^  -  3003^  +  3)  +  4B(aj^^  -  a^)pa^  -  1) 

+800^0^(703^  -  1)  +  2D(Oj^‘^  +  y 

2  2 

+8Eaj^Q;2(aj^  "  °2 

By  use  of  the  operator  equivalents  of  Table  IV-1  we  have  the  fourth  order 
part  of  the  spin  Hamiltonian 
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V^(S)  -  A  [3583“  -  308(8+1)83^  +  2583^  -  6S(8+l)+38^(8+l)^ 


3 


+b|(8^^  +  8_^)  1^783^  -  8(8+1)  -  5]  +[783^  -  S(8+l)-5] 

(8^  +  8,^1 

-iC  |(8^^  -  8,^)  1^783^  -  S(8+l)-5j  +  |^783^-8(8+1)-5] 

-  s.2)J 

+  D  +  8_'^J  -  IE  [^8^^  -  S_^]  , 


V-8 


where  8_|_  ■  8j^  +  83.  For  in  the  a-b  plane  a3  ■  0,  ■  cos  ^  and  03  ■  sin^, 

where  S  ia  the  angle  between  H  and  the  a  axis.  Therefore,  for  H  In  the 

7  Q  — .  '  O 

a-b  plane 

“  A^°/15  =  A  -  4B  cos  2j>  -  4c  sin  2^  +  2Dcos4^+2E  sin  4  (p  ,  V-9 


81milarly,  for  in  the  b-c  plane,  a3  >  cos  0,  ■  sin  9,  ■>  0  and  in 

this  plane 

-  A^°/15  -  I  (3A-4Bf2D)  +  (  |  A-2B-D)cos  29 
+(  +  7/2B  +  1/4D)  cos  49.  V-10 


The  best  fit  obtained  between  these  expressions  and  the  data  of  Fig.  13  was 
obtained  for 


A  -  -0.35 
B  -  -0.21 
-0.98 
D  ■  +0.53 
E  -  -0.59 


oe. 

oe. 

oe. 

oe. 

oe. 
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The  curves  calculated  from  these  values  are  plotted  In  Fig.  12  along 

with  the  data.  The  fit  Is  fair  or  good,  but  not  excellent.  It  should 

be  noted  that  the  data  may  contain  considerable  systematic  error,  since 

Che  points  were  calculated  from  smoothed  curves  drawn  through  the 

raw  data.  The  are  small  quantities,  on  Che  order  of  10  gauss,  obtained 

from  Che  subtraction  of  nearly  equal  resonance  fields  of  the  order  of 

10  gauss,  thus  magnifying  any  systematic  error. 

The  spin  Hamiltonian  coefficients  A  through  E  as  deduced  for 
3+ 

the  Gd  site  above  describe  a  potential  which  Is  very  far  from  cubic. 

For  a  cubic  potential  at  this  site  B,  C  and  E  would  be  zero,  whereas  In 
fact  C  Is  the  largest  coefficient  obtained,  and  E  Is  the  next  largest. 

Further,  for  a  cubic  site  D  =  -5A/2,  a  relationship  also  violated  in 

3+  3+ 

Che  Gd  data  fit.  In  short,  the  paramagnetic  resonance  results  on  Gd 

Indicate  Chat  the  rare  earth  Ion  site  in  the  orthoalumlnate  Is  so  severely 

distorted  from  the  cubic  prototype  that  no  visible  cubic Ity  remains. 

Calculations  of  crystal  field  effects  based  on  a  cubic  approximation 

would  appear  Invalid. 
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VI.  CRYSTAL  FIELD  CALCULATIONS. 

A.  The  Iron  Site. 

3+ 

The  dominant  crystal  field  effect  oh  an  S  state  Ion  like  Fe  Is 

usually  quadratic  In  the  spin  components.  This  quadratic  form  can  always 

be  diagonalized  by  a  proper  rotation  of  coordinates  and  the  direction  of 

the  principal  axes  and  their  magnitudes  specified.  The  resolution  of  the 

resonance  data  In  section  IV  has  given  us  the  orientation  and  shape  of 

3+ 

the  quadratic  energy  tensor  for  Fe  In  YAIO^. 

We  now  proceed  to  calculate  the  crystalline  electric  fields 
present  In  the  orthoferrite  on  a  point  Ion  model  and  compare  these  with 
the  data.  Since  the  mechanism  connecting  the  coefficients  In  the  spin 
Hamiltonian  and  the  crystal  fields  actually  present  Is  not  understood 
quantitatively  for  S  state  Ions,  we  cannot  pursue  a  truly  quantitative 
comparison.  We  can,  however,  calculate  the  quadratic  (In  coordinates) 
part  of  the  crystal  field  at  the  Iron  site  on  a  point  Ion  model  and 
compare  the  orientation  and  shape  of  this  quadratic  form  with  the 
experimentally  determined  spin  Hamiltonian. 

The  Iron  Ion  occupies  In  the  orthoferrite  a  position  of  sixfold 
oxygen  coordination,  as  illustrated  In  Fig.  13.  The  distortion  of  the 
coordination  octahedron  from  a  purely  cubic  configuration  Is  substantial. 
At  present  we  have  the  oxygen  parameters  only  for  GdFeO^;^  those  for 
YAIO^,  really  appropriate  for  this  calculation,  are  currently  being 
determined  by  Dr.  Ferdinand  Euler  at  Air  Force  Cambridge  Research  Center. 
Since  the  GdFeO^  parameters  are  the  only  ones  available  we  shall  therefore 


-  39  - 


proceed  using  them.  They  give  us  the  following  radius  vectors,  in  A,  from 
3+ 

the  Fe  ion  at  special  position  (1/2,  0,  0)  to  nearest  nel^bor  oxygens: 

r^  -  1.94  (-.139i'  -  .088^  +  .990k) 
rj  -  1.94  (+.1391  +  .088^ -  .990?) 
rjj  -  2.03  {-.nit  +  .620^  -  .187?) 

-  2.03  (+.772?  t  .620?  +  .187?) 
rj  -  1.94  (-.5771  -  .794? -  .196?) 
r^  -  1.94  (+.577?+  .794?+  .196?) 

where  ?,  ?  and  ?  are  unit  vectors  in  the  unit  cell  a,  b  and  c  directions 

respectively,  and  the  oxygens  are  Indexed  according  to  the  scheme 

Illustrated  in  Fig.  13.  The  potential  at  a  point  R6(^  within  the 

33 

coordination  cube  can  be  expressed  as 


00  m-4ti 


i*l  n"0  m»-n 


idiere  the  are  spherical  harmonics,  the  summation  over  1  is,  in  this 
case,  over  the  oxygen  neighbors,  and  the  >^^9^  ^  are  the  spherical 

coordinates  of  the  neighbors.  For  the  present  we  are  concentrating 
on  the  quadratic  portion  of  the  potential,  that  corresponding  to  n  ■  2. 
Performing  the  indicated  summations  for  n  -  2  one  obtains  the  potential 


VI- 1 


VI- 2 
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VI-3 


V(2)  -VpR^(  ~  jj.l 


i  2.6)Y2^ 


+  -  (8.4+i  2.6)Y,"V  -  (-4.3fl  4.2)  Y,'^^ 

/T  ^  /r  2 


Several  proportionality  constants  have  been  absorbed  Into  but  Its 
exact  definition  is  not  liiq>ortant  here  since  we  wish  to  examine  only  the 
location  and  relative  sizes  of  the  principal  axes  of  the  form.  To 
acconpllsh  the  dl agonal Izat ion  the  potential  was  first  put  into  the 
coordinate  form: 


V(2)  -  [-4.70  x^  -  0.40  y^  +  5.10  z^ 


-4.20  xy  +  2.60  yz  +  8.40  xz  J , 

The  transformation  T  was  then  sought  which  diagonalized  the  quadratic 
form;  this  same  transformation  yields  the  directional  cosines  of  the 
principal  axes  of  the  form.  The  desired  transformation  is: 


and  the  principal  values  are 

\(x')  -  -7.08 
\(y')  -  0.41 
\(z')  -  6.67 


VI-4 


VI-5 


VI- 6 
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Let  us  now  con^are  these  results  with  our  spin  Hamiltonian,  idiose  principal 
values  were  0^  ■  -1.45,  O2  •  0.27,  0^^  ■  1.18.  Superficially  these  are 
similar,  one  large  negative  value,  one  positive  value  almost  as  large,  and 
one  small  positive  value.  Correlation  ceases  at  this  point,  however.  The 
principal  axes  of  the  two  tensors  bear  little  relation  to  one  another. 

Only  one  close  correlation  exists.  In  that  the  axis  of  the  spin  Hamil¬ 
tonian  Is  close  to  the  y'  axis  of  the  point  Ion  potential,  but  this  Is  at 
best  an  unfortunate  correlation,  since  Is  associated  with  the  large 
negative  principal  value  idiereas  y’  Is  associated  with  the  small  positive 
principal  value. 

A  very  likely  source  of  this  lack  of  correlation  Is  the 

circumstance  that  calculations  based  on  the  oxygen  coordinates  for  the 

GdFeO^  are  used  to  compare  with  data  taken  on  YAIO^.  It  has  been  shown 

by  other  Investigators  that  the  quadratic  potential  Is  particularly 

35 

sensitive  to  crystalline  distortions.  Remal  and  deMars,  in  a  study 

3+ 

of  the  paramagnetic  resonance  spectrum  of  Gd  In  gallium  and  aluminum 
garnets,  showed  that  the  second  order  potential  coefficients  changed  by 
factors  of  from  two  to  four  in  going  from  YGaG  to  YAIG,  while  the  fourth 
order  potential  coefficients  were  changing  5  to  10  percent.  Only  further 
work,  using  the  positional  parameters  appropriate  for  YAIO^  can  resolve 
this  question. 

B.  The  Rare  Earth  Ion  Site. 

The  rare  earth  Ions  In  the  orthoferrite  occupy  an  unusual  site 
with  twelve  nearest  neighbor  oxygens,  which  may  be  visualized  as  occupying 
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the  midpoints  of  the  twelve  edges  of  a  cube  surrounding  the  rare  earth 
site.  The  elgjht  corners  of  the  cube  are  occupied  by  trlvalent  cations^ 

Fe  In  YFeO^j  or  A1  In  YAIO^.  In  the  orthorhombic ally  distorted  perov- 
sklte  the  cation  cube  remains  rigorously  cubic,  but  the  rare  earth 
Ion  Is  displaced  (In  the  a-b  plane)  from  the  center  of  this  cube  and  the 
oxygens  are  displaced  from  their  positions  at  the  midpoints  of  the  cube 
edges. 


We  have  not  calculated  the  crystal  fields  at  the  rare  earth  Ion 

for  the  distorted  perovsklte.  The  results  of  the  Immediately  previous 

section  Indicate  strongly  that  the  positional  parameters  obtained  from  GdFeO^ 

do  not  yield  useful  potentials  for  Interpreting  the  data  on  YAIO^.  We 

have,  however,  made  some  crystal  field  calculations  In  the  cubic 

c^proxlmatlon,  primarily  to  establish  an  expected  sequence  of  energy 
3+ 

levels  for  the  Yb  optical  studies.  The  potential  was  calculated  from 
expression  VI-2,  with  the  summation  In  1  now  running  over  the  twelve  oxygen 
neighbors.  The  resulting  potential  Is 

-4. 


JZOf^  +4 

4 


Y,  ■  ■(e(J>  )+Y, 
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j 

r-Vg°(e^)  +yi4  (0^)  +  Yg'^(e()))jj-  , 


VI-7 


where 

r  .  r  ^  ZeV 

C4  -  ^  9  )  5 

*  VI-8 

c  m  {  ^  ^1/2  ZeV 
^  13  ^  7 
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We  may  coflq>are  this  with  the  cubic  potentials  for  the  four,  six, 

36 

and  eightfold  coordinated  site,  as  calculated,  say,  by  Wolf,  et  al.  We 
find  that  the  coefficients  of  the  and  terms  are  of  the  same  relative 
sign,  as  In  the  case  for  the  octahedral  (six-fold  coordinated)  site,  but 
that  the  signs  of  both  the  Y^  coefficient  and  the  Y^°  coefficient  are 
opposite  for  the  12  fold  site  from  what  they  were  for  the  six  fold  site. 

We  therefore  expect  that  the  crystal  field  splitting  of  any  particular  J 
state  of  a  given  Ion  will  be  Inverted  In  the  twelvefold  case  as  conq>ared 
to  the  six  fold  case.  Specifically,  we  expect  that  the  excited  J  ■  5/2 
state  of  the  Yb  Ion  In  the  orthoferrite  will  have  the  quartet 
lowest  and  the  P  ^  doublet  higher;  and  that  for  the  J  ■  7/2  ground  state 
either  the  P ^  quartet  or  the  P  j  doublet  may  be  lowest. 

Referring  to  Fig.  14  we  observe  that  the  doublet  state  (  P^  7) 
lies  lower  than  the  quartet  (  ^  ?)  for  the  J  <■  5/2  state  of  Yb  ,  In 

direct  contradiction  to  the  predictions  of  the  preceding  paragraph. 

The  crystal  field  splitting  of  the  J  >  5/2  state  depends.  In 
the  cubic  approximation,  only  on  the  fourth  order  term  In  the  crystal 
field  potential.  The  P  ^  doublet  can  be  lowest  only  If  the  sign  of  the 
coefficient  of  Is  positive.  To  see  If  the  neglect  of  next- to- nearest, 
etc.  neighbors  could  be  affecting  the  sign  of  we  have  carried  this 
part  of  the  calculation  out  to  fourth  nearest  neighbors,  with  the  results 
shown  In  Table  VI-1.  The  coefficient  of  remains  negative  at  all  times. 
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Neighbor 

Radius 

Contribution 
to  potential 

Cumulative 
total  potential 

12  0^" 

a 

2_4 

■3.5^ 

-3.5  eV/a^ 

8  Al^ 

a  Un 

& 

2  4 

+3.39 

-.11  e^R^/a^ 

6 

a  (l 

a^ 

-1.86eV/a^ 

7  6. 

-1.97eV/a^ 

24  0^‘ 

a  JJ 

/e  2„4,  5 
-.45e  R  /a 

-2.42eV/a^ 

TABLE  VI- 1 


and  converges  on  a  negative  value  not  too  different  from  the  nearest 
neighbor  only  coefficient. 

We  cannot  get  the  correct  sequence  of  levels  for  the  5/2  state 

unless  we  abandon  the  cubic  approximation,  a  result  we  have  already 

3+ 

anticipated  from  the  Gd  spin  Hamiltonian. 
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VII.  ABSORPTION  SPECTRUM  OF  Yb^'*’ 

A.  Optical. 

3+ 

The  absorption  spectrum  of  Yb  In  the  near  Infrared  was  studied 
at  high  resolution  and  at  low  temperatures  (down  to  4.2°  K).  The  absorption 
was  studied  both  for  Yb  dllutely  substituted  Into  YAIO^  and  for  the 
fully  substituted  YbFeO^*  It  was  anticipated  that  an  understanding  of  the 

crystal  fields  present  would  be  obtained,  and  that  direct  observations  of 

3+  3+  37“39 

the  Yb  -  Fe  exchange  would  be  made. 

Unfortunately,  thou^  Interesting  and  suggestive  data  has  been 
obtained,  the  objectives  have  been  attained  In  only  a  limited  degree. 
Considerably  more  Information  remains  to  be  uncovered  In  future  spectroscopic 
work. 

Ytterbium,  though  electronically  the  sliq>lest  of  the  rare  earths 
substituents  In  the  orthoferrites,  does  not  always  exhibit  a  spectrum 
which  Is  easily  Interpreted.  First  of  all  there  is  not  enough  redundancy 
In  this  very  sparse  spectrum  to  allow  a  cross-check  of  assignments.  In 
addition,  the  vibrational-electronic  Interaction  Is  especially  strong  for 
ytterbium,  and  combination  (vlbronlc)  lines  often  rival  the  pure  electronic 
lines  In  Intensity.  This  makes  the  Identification  of  pure  electronic 
levels  difficult.  Fortunately  the  vlbronlc  lines  are  generally  somewhat 
more  diffuse  than  the  electronic  lines.  Furthermore  the  vlbronlc  lines 
seem  generally  to  be  suppressed  In  the  dllutely-doped  crystals.  Thus 
the  logical  procedure  Is  to  Identify  Insofar  as  possible  the  pure 
electronic  transitions  for  Yb  dllutely  Incorporated  In  YAIO^,  and  then  to 
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correlate  theae  results  with  the  spectrum  of  YbFeO^.  This  procedure  also 
allows  one  to  first  examine  crystal  field  effects  and  then  Introduce 
exchange. 

If  we  accept  as  pure  electronic  only  the  strongest  or  sharpest 
lines  In  the  spectrum  of  (Y  993^^  obtain  the  energy  level 

schemes  shown  In  Fig.  14.  The  spectra  are  very  similar  (except  for  scale) 
to  those  observed  for  Yb  In  the  garnet  structure^  and  one  Is  tempted  to 
make  the  same  level  assignments,  Indexed  on  a  cubic  scheme.  The  (nearly) 
cubic  Indexing  Is  also  suggested  by  doubling  of  those  states  assigned  to 
the  r’g  representation.  That  the  crystal  fields  have  changed  In  going 
from  the  alumlnate  to  the  ferrite  Is  Immediately  apparent;  the  upper 
state  (J  ■  5/2)  splitting  Increases  50  percent,  while  the  splitting 

of  the  J  ■  7/2  ground  state  decreases  25  percent. 

We  have  already  discussed  the  problem  of  reconciling  a  point 
Ion  potential  (cubic  approximation)  to  the  level  sequence  of  the  J  5/2 
state  (sec.  VI-B).  One  might  take  the  position  that  point  Ion  calculations 
are  meaningless,  and  attenq>t  an  enq>lrlcal  fit  of  and  to  the  data. 

This  approach  also  leads  to  dilemmas.  One  finds  that  to  fit  the  data 
for  both  excited  and  ground  states  It  Is  necessary  to  assume  (a)  that 
the  potential  Is  much  larger  than  the  V^;  (b)  that  Vg  changes  sign  In 
going  from  the  ortho- alumlnate  to  the  orthoferrite.  Both  of  these  assumptions 
are  unpalatable,  based  on  other  experience  In  rare  earth  spectra,  and  one 
Is  again  driven  toward  the  abandonment  of  a  cubic  approximation. 
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Very  Interesting,  though  as  yet  Inconel ete,  Information  on  the 

3+  3+ 

exchange  coupling  between  the  Yb  and  Its  Fe  neighbors  was  obtained 

from  the  YbFeO^  spectrum.  (In  the  case  of  the  garnets  a  very  full 

37-39 

characterisation  of  this  exchange  was  possible  ).  The  sharpest 

line  In  the  YAIO^  spectrum  (the  P ^ — ^^7  transition  In  Fig.  14)  becomes 

a  diffuse  doublet  In  YbFeO^  with  a  splitting  at  78°  K  of  about  6  cm 

This  splitting  Is  Insensitive  to  ten^erature  down  to  8°  K,  at  which 

temperature  the  moment  of  YbFeO^  Is  known^^  to  flop  from  the  £  axis 

to  the  a  axis.  Below  this  temperature  a  diffuse  four  line  pattern  Is 

observed  with  larger  splittings  between  adjacent  lines  (8,  11  and  8  cm  ^). 

No  detailed  Interpretation  of  the  spectrtjm  has  yet  been  made, 

but  a  number  of  observations  could  be  made.  The  splittings  are  surprisingly 

3+  3+ 

large  In  view  of  the  fact  that  the  Fe  spins  and  the  Yb  spins  are 
constrained  by  symmetry  considerations  to  be  nearly  perpendicular  to  each 
other.  The  Isotropic  part  of  the  rare  earth- Iron  exchange  will  therefore 
be  Ineffective  In  polarizing  the  rare  earth  Ions  or  In  producing  a  Zeeman 

3^ 

splitting  In  the  rare  earth  spectrum.  (In  GdFeO^,  where  both  Gd  and 
3+ 

Fe  are  nearly  Isotropic  S  states,  virtually  no  polarization  of  the 
Gd^  Is  produced^^  and  5"^  Is  the  same  as  for  YFeO^  down  to  4.2°  K) 

The  exchange  splitting  In  the  Yb  spectrum  must  be  accomplished  by  the 
antisymmetric  part  of  the  Yb-Fe  exchange,  and  this  antlsynraetrlc  part  appears 
to  be  nearly  as  large  as  the  symmetric  exchange  observed  In  Yb  Iron  garnet. 
This  Is  In  Itself  a  remarkable  conclusion. 
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One  would  expect  a  four  line  spectrum  from  the  transitions 
between  two  Zeeman  split  doublets,  so  Che  presence  of  only  two  lines 
above  8°  K  might  appear  puzzling.  However,  the  spectra  were  taken  on 
as-grown  platelets  with  Che  llg^C  propagating  along  Che  c  direction. 

Since  this  Is  also  Che  direction  (above  8°  K)  of  polarization  of  the 
Yb  Ions,  and  of  the  effective  field  polarizing  these  Ions,  It  Is 
probable  that  polarization  effects  elinninatcd  the  two  n  lines  and  left 
only  the  two  (5*  lines.  Below  8°  the  magnetic  polarization  is  In  a  plane  J.  to 
Che  c  axis  and  all  four  components  would  be  observed.  Which  splittings 
belong  to  which  doublet  could  be  deduced  from  the  tenperature  dependence 
of  absorption  Intensities.  Such  an  Intensity  study  and  further  polari¬ 
zation  studies  should  be  pursued  to  derive  the  latent  Information  In  the 
optical  spectrum. 

B.  Paramagnetic  Resonance. 

Crystals  of  YAIO^  doped  with  .5  percent  and  with  5  percent  Yb 

were  examined  at  35  kMc  and  at  9  kMc,  for  their  paramagnetic  resonance 

spectrum.  No  absorption  identifiable  as  arising  from  Yb  were  found 

at  magnetic  fields  up  to  20  kilogauss  or  at  temperatures  down  to 

1.4°  K,  or  for  a  wide  range  of  microwave  power  levels  or  field 

modulation  Intensities.  The  spectra  of  a  considerable  variety  of 

3+  3+ 

Impurities  were  found.  Including  Fe  ,  N1  and  Pc  In  an  unidentified 
40  3f 

valence  state.  Yb  was  known  to  be  present  In  the  samples  from 
Che  optical  spectra. 
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3f 

Die  paramagnetic  resonance  of  the  ground  state  doublet  of  Yb 

In  YAIO^  was  expected  to  yield  considerable  correlative  Information  for 

Interpreting  the  optical  absorption  spectrum  and  for  determining  the 

crystal  fields  present.  The  lack  of  a  spectrum  Is  quite  surprising  and 

strongly  suggests  that  the  groimd  state  Is  such  that  ■  1  transitions 

between  the  two  partners  of  the  lowest  doublet  are  forbidden.  Such  can 

happen^  for  Instance,  for  an  axial  field  which  might  put  a  J  >  +  7/2 

state  lowest.  It  cannot  happen  for  a  cubic  field.  The  negative  results 

of  the  resonsnce  search  tend  to  Indicate  that  the  Pj  doublet  does 

not  lie  lowest  In  the  orthoferrite,  since  the  paramagnetic  characteristics 

41  42 

of  this  state  are  well  known.  * 

It  Is  also  possible  that  variations  In  the  local  geometry,  due 
to  some  defect  common  In  the  orthoferrite  structure,  broadens  the 
paramagnetic  resonance  so  much  as  to  make  It  unobservable.  There  Is 
some  evidence  (satellite  lines)  In  the  optical  spectra  of  considerable 
site  Irregularity.  It  Is  hard  to  estimate  how  seriously  to  take  the 
line -broadening  explanation  of  the  negative  results. 
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VII  SUMMARY  AND  CONCLUSIONS. 

A  number  of  pieces  will  have  to  be  fit  together  before  one 
can  understand  the  macroscopic  magnetic  properties  of  the  orthoferrites 
from  an  atomic  point  of  view.  These  pieces  of  Information  are  comprised 
of  knowledge  of  the  low-lying  levels  of  the  magnetically  active  lons^ 
the  relationship  of  these  levels  to  the  Ion  and  Its  crystalline  environ¬ 
ment^  and  the  Interaction  between  Ions. 

We  have  In  the  course  of  this  contract  work  assembled  a  number 
of  these  pieces,  but  far  too  few  to  complete  the  picture;  and  It  would 
seem  that  we  have  also  collected  some  pieces  of  Information  belonging  to  a 
different  puzzle.  In  particular  It  has  become  apparent  that  much  of 
the  detailed  Information  collected  In  the  ortho-aluminate  cannot  be 
applied  without  considerable  interpretation  to  the  orthoferrites. 

We  have  accomplished  the  following; 

(1)  Characterized  the  spin  configurations  allowed  In  the 
orthoferrites  from  magnetic  group  theory. 

(2)  Derived  the  antiferromagnetic  resonance  frequencies, 
static,  and  rf  susceptibilities  of  yttrium  orthoferrite. 

(3)  Measured  the  paramagnetic  resonance  spectra  of  Fe  and 
Gd^  in  YAlOg. 

(4)  Developed  a  powerful  method  for  extracting  the  spin 
Hamiltonian  for  S  state  Ions  on  sites  of  low  symmetry. 

(5)  Solved  for  the  spin  Hamiltonian  of  Fe  and  Gd  ,  finding 
that  the  Fe  is  nearly  cubic,  but  that  the  rare  earth 
site  Is  so  distorted  that  no  cublclty  remains. 
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3+ 

(6)  Concluded  from  the  Fe  spectrum  that  single  Ion 
anisotropy  effects  would  not  be  large  enough  to  cause  the 
observed  canting  In  YFeO^. 

3+ 

(7)  Compared  the  spin  Hamiltonian  of  Fe  In  YAIO^  with 
crystal  field  calculations  based,  on  GdFeO^  oxygen 
parameters,  and  found  little  correlation. 

3+ 

(8)  Found  on  the  basis  of  optical  spectroscopy  (of  Yb  ) 
that  the  crystal  fields  In  YAlO^  are  substantially 
different  from  those  In  YbFeO^. 

(9)  Found  from  optical  observations  of  exchange  splittings 
In  YbFeO^  that  a  strong  antisymmetric  exchange  (S^  x  S2) 
exists  between  the  Yb  and  Fe  magnetizations. 

(10)  Indicated  a  number  of  directions  for  further  work  which 
may  clear  up  the  questions  and  contradictions  raised 
by  the  Investigations  to  date. 
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Unit  cell  of  the  orthorhombic  orthoferrite,  showing  oxygen 
and  Iron  sites. 


o  RARE  EARTH 


FIGURE  3 


Unit  cell  of  the  orthorhombic  orthoferrite,  showing  and 
Indexing  the  rare  earth  Ion  sites. 
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FIGURE  4 

3+ 

Paramagnetic  resonance  absorption  spectrum  of  Fe  In 
YAlO.,  with  H  In  the  a-b  plane. 
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IN  DEGREES  BETWEEN  H  AND  b  AXIS 

FIGURE  5 

3+ 

resonance  absorption  spectrum  of  Fe  In 
in  the  b-c  plane. 


IN  KILOGAUSS 


FIGURE  7 

Pl<^  of  Hamiltonian  parameter  A,”  aa  a  function  of  angle  for 
Fe^  In  YAIO3.  The  angle  la  meaaured  from  the  b,  b,  and 
c  axea  In  the  a-b.  b-c  and  apecial  planea,  reapMtTvely. 
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ANGLE  IN  DEGREES  BETWEEN  H  AND  [a]  AXIS 


FIGURE  9 

3 

Paramagnetic  resonance  absorption  spectrum  of  Gd 
YAIO3,  with  In  the  a-b  plane. 
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FIGURE  12 

Plot  of  Hamiltonian  parameter  as  a  function  of  angle 
for  Gd^'*'  in  YAIO3.  The  angle  is  measured  from  the  a  and 
c  axes  in  the  a-b  and  b-c  plmes,  respectively. 
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OXYGEN 


#  IRON 


FIGURE  13 


Coordination  polyhedron  for  Fe  In  YAl03y  showing  the 
positions  of  the  six  nearest  neighbor  oxygens. 
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